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Abstract 

Electronic states in nanographite ribbons with zigzag edges and zigzag carbon nanotubes are 
studied using the extended Hubbard model with nearest neighbor Coulomb interactions. The 
nearest Coulomb interactions stabilize electronic states with the opposite electric charges sepa- 
rated and localized along both edges. Such states are analogous as nanocondensers. Therefore, 
electric capacitance, defined using a relation of polarizability, is calculated to examine nano 
functionalities. We find that the behavior of the capacitance is widely different depending on 
whether the system is in the magnetic or charge polarized phases. In the magnetic phase, the 
capacitance is dominated by the presence of the edge states while the ribbon width is small. 
As the ribbon becomes wider, the capacitance remains with large magnitudes as the system 
develops into metallic zigzag nanotubes. It is proportional to the inverse of the width, when the 
system corresponds to the semiconducting nanotubes and the system is in the charge polarized 
phase also. The latter behavior could be understood by the presence of an energy gap for 
charge excitations. 
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1. Introduction 



Nano carbon materials have been attracting much attention both in the fundamental science 
and in the interests of application to nanotechnology devices. 1,2 > Their physical and chemical 
natures change variously depending on geometries. 1-3 ) In carbon nanotubes, diameters and 
chiral arrangements of hexagonal pattern on tubules decide whether they are metallic or not. 1 ' 2 ) 
In addition, carbon systems with open boundaries show unusual features. 

In nanographites, the peripheral edges strongly affect the electronic states, 3 ) and the nanographite 
ribbon with zigzag edges has the nonbonding molecular orbitals localized mainly along the 
zigzag edges. The partly flat electronic bands give rise to a sharp peak in the density of states. 
The on-site Coulomb interaction by the Hubbard model has suggested to induce the spin po- 
larized (SP) state, and the most spin density concentrates along the zigzag edges forming the 
ferrimagnetic orders due to the edge states. As to the effects of the Coulomb interactions, 
another qualitatively important problem has come out. This is by no means the competition 
between the spin and charge orderings in the nanographite system due to the on-site and near- 
est neighbor Coulomb interactions. In fact, for two dimensional (2D) boundless graphite, the 
possibility of the stable charge density wave (CDW) state compared with the spin density wave 
(SDW) state was reported. 4,5 ) 

Recently, 6 ^ we have demonstrated that the nearest neighbor Coulomb interaction stabilizes 
a novel charge polarized (CP) state with a finite electric dipole moment in zigzag ribbons, and 
that it competes with the SP state. Such the CP state is caused by the interplay between the 
nearest neighbor Coulomb interaction and the edge states in nanographene ribbons with zigzag 
edges using the extended Hubbard model. This study has been treated rather formally by 
changing the interaction parameters in the wide ranges. Though it has been discussed that the 
transverse electric field might induce the first order phase transition from the SP state to the 
CP state, 6 ) we need further study in order to reveal what roles such the novel SP and CP states 
play in actual physical quantities measurable in experiments. Such theoretical studies should 
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be useful when materials specific properties due to edge localized states can be expected. 

In this paper, electric capacitance is calculated extensively in order to examine possible 
functionalities as "nano-size condensers" of nanographite materials. The positive and negative 
charge accumulations along the edges of the both sides have a pattern analogous to a condenser, 
if each set of edge site atoms is regarded as one electrode. One of them is the positive electrode, 
and the other is the negative one. The spacer between two electrodes is the inner part of the 
zigzag nanoribbon. The capacitance is calculated with the relation of the charge polarizability. 
This definition is the usual capacitance for the SP state because there is not charge polarization 
without an electric field. However, the definition becomes different for the CP state, because 
the self-polarized charge is present at the edge sites. In this case, the calculated capacitance is 
a "differential capacitance" as defined by the relation of polarizability. 

The calculated electric capacitance will be discussed with scopes of nano functionalities. 
In the magnetic phase, the capacitance is dominated by the presence of the edge states while 
the ribbon width is small. As the ribbon becomes wider, the capacitance remains with large 
magnitudes as the system develops into metallic zigzag nanotubes, while it is proportional to 
the inverse of the width when the system corresponds to the semiconducting nanotubes. In the 
charge polarized phase, the capacitance is inversely proportional to the graphite width. 

This paper is organized as follows. In section 2, we explain the extended Hubbard model. 
In section 3, we give a brief review of the SP and CP states in order to look at electronic 
properties. In section 4, calculated electric capacitance is reported and discussion is given. The 
paper is closed with a short summary in section 5. 

2. Model and method 

Figure 1 illustrates the geometry of a zigzag ribbon with an inversion symmetry. Here, N 
and L are the width and length of the ribbon, respectively. The dashed rectangle denotes a 
unit cell. The periodic boundary condition is set along the y axis of Fig. 1. Since the zigzag 
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ribbon is bipartite, A and B sites are assigned by filled and open circles, respectively. All the 
twofold coordinated sites in the lower and upper zigzag edges belong to the A and B sublattices, 
respectively, at which edge states are mainly localized. 

We consider the 7r-electron system on the zigzag ribbon using the extended Hubbard Hamil- 
tonian with the on-site U and nearest-neighbor V Coulomb interactions. The model is as 
follows: 

H = -t £ (cU-,,, + h.c.) + £/£(4 T c, T - ^-)(clc hl - ^) 

(i,j),o- i 
+ V Y,(Y, C laCi,a-n cl )(^2 C ],r C j,r-n cl ), (1) 

where t (> 0) is the hopping integral between the nearest neighbor ith and jth sites; the sum 
with is taken for all the pairs of the nearest neighbor sites; c i>a annihilates an electron of 
spin a at the ith site; n e \ is the average electron density of the system. 

We adopt the standard Unrestricted Hartree-Fock approximation to this model: 3,6 ' 7 ) 

4t c *>t c L c a <4,T c ^T> c I,i c ^i + c ha,i(c{i c id) - ( c h c ^)( c {i c i,i)- ( 2 ) 
3. Magnetic and charge ordered states 

In this section, we briefly review the competition between the SP and CP states, in order to 
give an introduction to the materials and their electronic properties. 

Figure 2 (a) shows the charge density distribution in the CP state with iV x L = 4 x 20, 
where U = 0.3t and V = OAt, t being the hopping integral between neighboring carbons. For 
comparison, Fig. 2 (b) shows the spin density profile of the SP state for U — l.Oi and V — 0. 
The CP state has no spin density at every site, while the SP state has no charge density at 
every site. In the CP state, the upper zigzag edge is charged positively, while the lower zigzag 
edge is charged negatively. The distribution pattern of the charge density in the CP state is 
just like that of the spin density in the SP state. The charge polarization in the CP state 
can be explained by the interplay between the edge states and the Fermi-instability of the flat 
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bands due to the interaction V. The CP state has a finite electric dipole moment pointing from 
the upper edge to the lower edge along the x-axis of Fig. I. It is formed in the homogeneous 
system and breaks the inversion symmetry. 

The relative stabilities between the CP and SP states have been investigated and summa- 
rized in the phase diagram. 6 ) While U is fairly larger than V , the SP state is more stable. As 
V becomes stronger, the system exhibits the first order phase transition from the SP state to 
the CP state. The charge order of Fig. 2 (a) occurs in the CP phase, while the spin order of 
Fig. 2 (b) in the SP phase. 

4. Electric capacitance 

In order to examine nano functionalities as"nano-size condensers", electric capacitance of the 
nanographite ribbons is calculated. The idea of the calculation is shown schematically in Fig. 
3. We assume that the two sets of carbon atoms at the zigzag edges are regarded as positive and 
negative electrodes, respectively. The absolute value of the net variation of the accumulated 
charge is divided by the strength of the small applied voltage, and the capacitance is obtained. 
In Fig. 3 (b), the quantity Qo is the net charge over the L/2 carbon atoms at the zigzag edge 
sites, when the static electric field is absent. The part dQ is the change of the net charge with 
respect to the small field which is parallel to the z-axis as shown in Fig. 3 (a). The capacitance 
C is calculated using the relation of polarizability dQ = CdV, where dV is the change of voltage 
due to the static electric field. The similar method has been used for the calculation of the 
capacitance using the first principle method. 8 ) We take dV = O.Olt and t = 2.5 eV. We also 
assume that the bond length between carbons is 1.45 A. 

Figure 4 shows the calculated results for the system in the SP state. The actual magnitudes 
in the logarithmic scale are plotted in Fig. 4 (a) against the ribbon width in the scale of A, and 
their inverse values are shown in (b). When the width iV of the system increases, it develops 
into the (L/2, 0) zigzag nanotube due to the periodic boundary condition along the y direction 
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of Fig. I. The nanotube is metallic when L/2 is a multiple of three, and is semiconducting for 
others 2 ^ in the noninteracting model. While the ribbon width is small (approximately < 10 A), 
the capacitance would be dominated by the presence of the edge states. Even though the ribbon 
width becomes larger, the capacitance remains with large magnitudes for L = 18, as shown by 
squares in Fig. 4 (a). This would be related with the resulting metallic properties of the long 
enough (9, 0) nanotube. On the other hand, the capacitance for L = 20 is proportional to the 
inverse width for larger widths, as shown in Fig. 4 (b). The presence of the semiconducting 
gap in the (10, 0) nanotube might result in the distinctive difference. 

Figure 5 summaries the calculated capacitance for the system in the CP state. The two 
sets of L = 18 and 20 are shown, too. We find that the capacitance is inversely proportional 
to the graphite width [Fig. 5 (b)]. This behavior does not depend on whether the long enough 
system is the metallic (9,0) or semiconducting (10,0) nanotube. As the system has the charge 
orders, electrons would become less mobile than those in the system without ordered states. 
The behavior could be understood by the presence of an energy gap for charge excitations. 

We look at interaction parameter dependences of the calculated capacitance. Figures 6 
and 7 show the results for the ribbon length L — 18 and 20, respectively. When the ribbon 
width becomes wider, the system develops into (9,0) and (10,0) zigzag nanotubes. The former 
is metallic and the latter is semiconducting. The Coulomb interaction strengths are (U, V) = 
(1.0*, 0) and (2.0t, 0) for the SP state, while (U, V) = (0, 0.5*) and (0, 0.6*) for the CP state. As 
the interaction strength U increases, the spin orderings become stronger. This gives rise to a 
stiff response to the applied field, thus the magnitude of C decreases in the SP state. Similarly, 
the magnitude of C decreases in the CP state, as the charge orderings become stronger with 
the increase of the parameter V. 

The intrinsic capacitance of the single wall carbon nanotube can be estimated as follows. 
We consider a nanotube with finite length which is present in a quantum box of the length L y . 
The wavenumber is quantized with the interval Ak = 2n/L y . Addition energy of an electron 
to the nanotube in the quantum box can be equated with an energy as a dielectric system. 
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The addition energy can be estimated to be AE = (h/2ir)vFAk ■ (1/2) • (1/2), where vf is the 
Fermi velocity, the first factor (1/2) comes from the degeneracies of spin, and the second one 
(1/2) is due to the degeneracy of electronic states near the Fermi energy. Here, we assume 
that the nanotube interacts with circumstances, and the degeneracies of spin and electronic 
states are lifted slightly. Therefore, AE is the overall average level spacing including spin and 
number of electronic bands. By setting AE = e 2 /2CqL y , we obtain the quantum capacitance 
per unit length: Cq = 2e 2 /hvF- For carbon nanotubes 

Vf = 8 x 10 5 (m/s), 9 ) so that C Q = 100 
(aF//x m) = 10~ 20 (F/A). In fact, the capacitance per unit length has been measured to be 190 
(aF//i m) 10 ) for example, and the above rough estimation explains the experimental magnitude 
fairly well. Such the order of magnitudes of the estimation and experimental value also agrees 
with that of the capacitance obtained in the calculation within the difference of a few order of 
magnitudes. Even though the detailed definitions of the capacitances are different mutually, 
the resulting order of magnitudes should reflect quantum characters of electronic systems. 

We can give a reason of the obtained magnitudes of the capacitance by analogy with a 
classical nano-size condenser, too. The calculated order of magnitudes of the capacitance 
around 10~ 20 ~ 10~ 19 F might be estimated as follows. When we consider a classical parallel- 
electrode condenser with the area of the electrode S, the capacitance is given by C = e(S/d), 
where e is a dielectric constant of the spacer, and d is the distance between the electrodes. We 
shall assume that S for the nanographite ribbon is estimated as (the effective thickness of the 
graphite plane ~ the interlayer distance of the stacked graphite ~ 3.5 A) x (the length of the 
ribbon). We obtain the magnitude C ~ 3 x 10 -20 F for iV = 2, L = 20, and e of the empty 
space. Though uncertainties of a few factors in e will be present as a material constant, the 
similar order of magnitudes would be obtained. Therefore, the calculated order of magnitudes 
seems reasonable for nano-size systems. 

The inversely proportional behavior against the width might indicate that electric field 
between the two electrodes is confined within the graphite plane, and does not come out to the 
outside. If the electric field comes out and the each electrode can be regarded as a metallic 
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wire, the capacitance would show a logarithmic dependence as a function of the ribbon width. 

5. Summary 

In nanographite ribbons with zigzag edges, 6 ) the charge-polarized state is present. This state has 
a finite electric dipole moment caused by the interplay between the nearest neighbor Coulomb 
interaction and the edge states, and competes with the spin-polarized state caused by the 
on-site Coulomb interactions. 

In this paper, electric capacitance has been calculated in order to test the nano functional- 
ities. In the magnetic phase, the capacitance is dominated by the presence of the edge states 
while the ribbon width is small. As the ribbon becomes wider, the capacitance remains with 
large magnitudes as the system develops into metallic zigzag nanotubes, while it is proportional 
to the inverse of the width when the system corresponds to the semiconducting nanotubes. In 
the charge polarized phase, the capacitance is inversely proportional to the graphite width. 
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Figure Captions 



Fig. 1. Schematic structure of the nanographite ribbon with zigzag edges (zigzag ribbon). The 
filled and open circles are A and B sites, respectively. The dashed rectangle denotes a unit cell. 

Fig. 2. (a) Charge density distribution of the charge-polarized (CP) state for U = 0.3t and 
V = OAt, and (b) the z component of spin density distribution of the spin-polarized (SP) state 
for U = t and V = 0, on a zigzag ribbon with 4 x 20 sites. The open and filled circles show 
positive and negative densities. Their radii are proportional to the magnitudes of the charge 
or spin densities: the maximum is 0.18 in (a) and 0.12 in (b). 

Fig. 3. (a) Periodic nanoribbon. Periodic boundary condition is applied to the one- dimensional 
direction of the nanoribbon. The nanoribbon strip is perpendicular to the x-y plane. Small 
static electric field is parallel to the z-axis. (b) Analogy with a parallel electrode condenser. 
The quantity Q is the net charge over the L/2 carbon atoms at the zigzag edge sites, when 
the static electric field is absent. The part dQ is the change of the net charge with respect to 
the small field. The capacitance C is obtained using the relation of polarizability dQ = CdV, 
where dV is the change of voltage due to the static electric field. 

Fig. 4. The electric capacitance calculated for the SP state at (U,V) = (1.0t,0). Two sets 
of the ribbon lengths L — 18 (squares) and 20 (circles) are considered. The magnitude of the 
capacitance (a) and its inverse (b) are plotted against the ribbon width in the scale of A. 

Fig. 5. The electric capacitance calculated for the CP state at (U,V) = (0,0. 5t). Two sets 
of the ribbon lengths L — 18 (squares) and 20 (circles) are considered. The magnitude of the 
capacitance (a) and its inverse (b) are plotted against the ribbon width in the scale of A. 
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Fig. 6. The electric capacitance calculated for the ribbon length L = 18. The Coulomb 
interaction strengths are (U, V) = (l.Ot, 0) and (2.0t, 0) for the SP state, while (U, V) = (0, 0.5*) 
and (0, 0.6t) for the CP state. The magnitude of the capacitance (a) and its inverse (b) are 
plotted against the ribbon width in the scale of A. 

Fig. 7. The electric capacitance calculated for the ribbon length L = 20. The Coulomb 
interaction strengths are (U, V) = (l.Ot, 0) and (2.0t, 0) for the SP state, while (U, V) = (0, 0.5*) 
and (0, 0.6*) for the CP state. The magnitude of the capacitance (a) and its inverse (b) are 
plotted against the ribbon width in the scale of A. 
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